A dynamic in vitro developing testis model reflects structures and functions of testicular development in vivo.

To better define appropriate applications of our 3-dimensional testicular co-culture as a model for reproductive toxicology, we evaluated the ability of the model to capture structural and functional elements that can be targeted by reproductive toxicants. Testicular co-cultures were prepared from postnatal day 5 male rats and cultured with a Matrigel overlay. Following a 2-day acclimation period, we characterized functional pathway dynamics by evaluating morphology, protein expression, testosterone concentrations, and global gene expression at a range of timepoints from experimental days 0-21. Western blotting confirmed expression of Sertoli cell, Leydig cell, and spermatogonial cell-specific protein markers. Testosterone detected in cell culture media indicates active testosterone production. Quantitative pathway analysis identified Gene Ontology biological processes enriched among genes significantly changing over the course of 21 days. Processes enriched among genes significantly increasing through time include general developmental processes (morphogenesis, tissue remodeling, etc.), steroid regulation, Sertoli cell development, immune response, and stress and apoptosis. Processes enriched among genes significantly decreasing over time include several related to male reproductive development (seminiferous tubule development, male gonad development, Leydig cell differentiation, Sertoli cell differentiation), all of which appear to peak in expression between days 1 and 5 before decreasing at later timepoints. This analysis provides a temporal roadmap for specific biological process of interest for reproductive toxicology in the model and anchors the model to sensitive phases of in vivo development, helping to define the relevance of the model for in vivo processes.

High-throughput screening tools facilitate calculation of a combined exposure-bioactivity index for chemicals with endocrine activity.

High-throughput and computational tools provide a new opportunity to calculate combined bioactivity of exposure to diverse chemicals acting through a common mechanism. We used high throughput in vitro bioactivity data and exposure predictions from the U.S. EPA's Toxicity and Exposure Forecaster (ToxCast and ExpoCast) to estimate combined estrogen receptor (ER) agonist activity of non-pharmaceutical chemical exposures for the general U.S. population. High-throughput toxicokinetic (HTTK) data provide conversion factors that relate bioactive concentrations measured in vitro (µM), to predicted population geometric mean exposure rates (mg/kg/day). These data were available for 22 chemicals with ER agonist activity and were estimated for other ER bioactive chemicals based on the geometric mean of HTTK values across chemicals. For each chemical, ER bioactivity across ToxCast assays was compared to predicted population geometric mean exposure at different levels of in vitro potency and model certainty. Dose additivity was assumed in calculating a Combined Exposure-Bioactivity Index (CEBI), the sum of exposure/bioactivity ratios. Combined estrogen bioactivity was also calculated in terms of the percent maximum bioactivity of chemical mixtures in human plasma using a concentration-addition model. Estimated CEBIs vary greatly depending on assumptions used for exposure and bioactivity. In general, CEBI values were <1 when using median of the estimated general population chemical intake rates, while CEBI were ≥1 when using the upper 95th confidence bound for those same intake rates for all chemicals. Concentration-addition model predictions of mixture bioactivity yield comparable results. Based on current in vitro bioactivity data, HTTK methods, and exposure models, combined exposure scenarios sufficient to influence estrogen bioactivity in the general population cannot be ruled out. Future improvements in screening methods and computational models could reduce uncertainty and better inform the potential combined effects of estrogenic chemicals.

Anchoring a dynamic in vitro model of human neuronal differentiation to key processes of early brain development in vivo.

We characterize temporal pathway dynamics of differentiation in an in vitro neurotoxicity model with the aim of informing design and interpretation of toxicological assays. Human neural progenitor cells (hNPCs) were cultured in differentiation conditions up to 21 days. Genes significantly changed through time were identified and grouped according to temporal dynamics. Quantitative pathway analysis identified gene ontology (GO) terms enriched among significantly changed genes and provided a temporal roadmap of pathway trends in vitro. Gene expression in hNPCs was compared with publicly available gene expression data from developing human brain tissue in vivo. Quantitative pathway analysis of significantly changed genes and targeted analysis of specific pathways of interest identified concordance between in vivo and in vitro expression associated with proliferation, migration, differentiation, synapse formation, and neurotransmission. Our analysis anchors gene expression patterns in vitro to sensitive windows of in vivo development, helping to define appropriate applications of the model.

High-throughput in Vitro Data To Inform Prioritization of Ambient Water Monitoring and Testing for Endocrine Active Chemicals.

The presence of industrial chemicals, consumer product chemicals, and pharmaceuticals is well documented in waters in the U.S. and globally. Most of these chemicals lack health-protective guidelines and many have been shown to have endocrine bioactivity. There is currently no systematic or national prioritization for monitoring waters for chemicals with endocrine disrupting activity. We propose ambient water bioactivity concentrations (AWBCs) generated from high throughput data as a health-based screen for endocrine bioactivity of chemicals in water. The U.S. EPA ToxCast program has screened over 1800 chemicals for estrogen receptor (ER) and androgen receptor (AR) pathway bioactivity. AWBCs are calculated for 110 ER and 212 AR bioactive chemicals using high throughput ToxCast data from in vitro screening assays and predictive pathway models, high-throughput toxicokinetic data, and data-driven assumptions about consumption of water. Chemical-specific AWBCs are compared with measured water concentrations in data sets from the greater Denver area, Minnesota lakes, and Oregon waters, demonstrating a framework for identifying endocrine bioactive chemicals. This approach can be used to screen potential cumulative endocrine activity in drinking water and to inform prioritization of future monitoring, chemical testing and pollution prevention efforts.

Identifying reference chemicals for thyroid bioactivity screening.

Reference chemicals were selected based on thyroid bioactivity in 'Tier 1' screening assays used by the U.S. EPA's Endocrine Disruptor Screening Program. Active reference chemicals had significant effects on thyroid-responsive endpoints in the amphibian metamorphosis assay, and the male and female pubertal rat assays. In the absence of thyroid weight or histopathological effects, additional published studies providing mechanistic data on thyroid activity were required for active chemicals. Inactive reference chemicals had no significant effects on thyroid-responsive endpoints in Tier 1 assays, or in amphibian or rodent studies from several online databases. The 34 reference chemicals (29 active and five inactive) will be useful for performance-based validation of alternative, high throughput screening assays for thyroid bioactivity.

The presence of macrophages and inflammatory responses in an in vitro testicular co-culture model of male reproductive development enhance relevance to in vivo conditions.

Our 3-dimensional testis co-culture system (3D-TCS) represents a promising model of male reproductive toxicity which captures sensitive processes of male reproductive development and contains the main testes cell types (germ, Leydig and Sertoli cells). Macrophages are another cell type important for testicular function and help to modulate immuno-endocrine processes during testes development. Chemicals such as phthalate esters (PE's) affect macrophage function and testosterone production in the testes in vivo. The aim of this study was to determine whether macrophages were present in the 3D-TCS and investigate responses in our model that may be related to immuno-endocrine functions. We observed consistent expression of the resident macrophage marker ED2 as well as increases in inflammatory cytokines produced by macrophages and testes cells (IL-6, TNF-α and KC/GRO) after exposure to toxic PE's. Pathway analysis of gene expression changes after exposure to PE's showed that IL-6 and TNF-α signaling pathways were enriched after treatment with reproductively toxic, but not non-reproductively toxic phthalates. These results indicate that macrophages and inflammatory processes are captured in the 3D-TCS and that these processes are impacted by exposure to reproductive toxicants. These processes represent a major mode of action for in vivo testis toxicity for a variety of compounds and our novel in vitro model is able to capture toxicant perturbation of immune function.

Differential epigenetic effects of chlorpyrifos and arsenic in proliferating and differentiating human neural progenitor cells.

Understanding the underlying temporal and mechanistic responses to neurotoxicant exposures during sensitive periods of neuronal development are critical for assessing the impact of these exposures on developmental processes. To investigate the importance of timing of neurotoxicant exposure for perturbation of epigenetic regulation, we exposed human neuronal progenitor cells (hNPCs) to chlorpyrifos (CP) and sodium arsenite (As; positive control) during proliferation and differentiation. CP or As treatment effects on hNPCs morphology, cell viability, and changes in protein expression levels of neural differentiation and cell stress markers, and histone H3 modifications were examined. Cell viability, proliferation/differentiation status, and epigenetic results suggest that hNPCs cultures respond to CP and As treatment with different degrees of sensitivity. Histone modifications, as measured by changes in histone H3 phosphorylation, acetylation and methylation, varied for each toxicant and growth condition, suggesting that differentiation status can influence the epigenetic effects of CP and As exposures.

Phthalate metabolism and kinetics in an in vitro model of testis development.

We have developed an in vitro model of testis development (3D-TCS) using rat testicular cells overlaid with extracellular matrix. One barrier preventing utilization of in vitro models in toxicity testing is the absence of metabolic capability. Another challenge is lack of kinetic data for compounds in vitro. We characterized metabolic capabilities and investigated the kinetics of phthalate male reproductive toxicants in the 3D-TCS. Cells were treated with three phthalate diesters for 2, 8 and 24 h. Parent compounds and metabolites were measured in cell culture media and cell lysate via mass spectrometry. Levels of monoester metabolites were used as an indication of metabolism of phthalates via lipase activity. Metabolites were detected in all treated cell media and cell lysate samples, with levels ranging from <0.5-14.7% of initial mass of parent compound. Phthalates partitioned between media and lysate in a manner consistent with each compound's degree of lipophilicity. UDGPT activity was detected in DBP and DEP treated samples. 3D-TCS microarray data indicated gene expression for lipases and CYPP450s. Results indicate that the 3D-TCS is a metabolically active co-culture and that physiochemical properties can provide information about the kinetics of compounds in the 3D-TCS, improving our ability to interpret results from the model.

Stage-specific signaling pathways during murine testis development and spermatogenesis: A pathway-based analysis to quantify developmental dynamics.

Shifting the field of developmental toxicology toward evaluation of pathway perturbation requires a quantitative definition of normal developmental dynamics. This project examined a publicly available dataset to quantify pathway dynamics during testicular development and spermatogenesis and anchor toxicant-perturbed pathways within the context of normal development. Genes significantly changed throughout testis development in mice were clustered by their direction of change using K-means clustering. Gene Ontology terms enriched among each cluster were identified using MAPPfinder. Temporal pathway dynamics of enriched terms were quantified based on average expression intensity for all genes associated with a given term. This analysis captured processes that drive development, including the peak in steroidogenesis known to occur around gestational day 16.5 and the increase in meiosis and spermatogenesis-related pathways during the first wave of spermatogenesis. Our analysis quantifies dynamics of pathways vulnerable to toxicants and provides a framework for quantifying perturbation of these pathways.

Preparation of rodent testis co-cultures.

Male reproductive development is a complex process that is sensitive to disruption by a range of toxicants. There is a great need for in vitro models that can evaluate potential male reproductive toxicants. The current unit presents a protocol for preparation of a three-dimensional in vitro model of male reproductive development that reduces the number of animals required for evaluation of toxicants. A Matrigel overlay provides a three-dimensional extracellular matrix that improves cell attachment, viability, and communication, and makes the model more reflective of in vivo environments.